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Observing the onset of outflow collimation in a massive protostar 
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The current paradigm of star formation through accretion disks, and magnetohydro- 
dynamically driven gas ejections, predicts the development of collimated outflows, rather 
than expansion without any preferential direction. We present radio continuum observa¬ 
tions of the massive protostar W75N(B)-VLA 2, showing that it is a thermal, collimated 
ionized wind and that it has evolved in 18 years from a compact source into an elongated 
one. This is consistent with the evolution of the associated expanding water-vapor maser 
shell, which changed from a nearly circular morphology, tracing an almost isotropic out¬ 
flow, to an elliptical one outlining collimated motions. We model this behavior in terms of 
an episodic, short-lived, originally isotropic, ionized wind whose morphology evolves as it 
moves within a toroidal density stratification. 


Water-vapor masers at 22 GHz are commonly found in star-forming regions, arising in 
the shocked regions created by powerful outflows from protostars in their earliest phases of 
evolution ( 1 ). Observations of these masers with very long baseline interferometry (VLBI) 
indicate that at the early life of massive stars there may exist episodic, short-lived (tens of 
years) events associated with very poorly collimated outflows ( 2 , 3 , 4 , 5 ). These results are 
surprising because, according to the core-accretion model for the formation of massive stars 
(>8 M 0 ), which is a scaled-up version of low-mass star formation, collimated outflows are 
already expected at their very early phases ( 6 , 7 , 8 ). 
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A unique case of a short-lived, poorly collimated outflow is the one found in the high- 
mass star-forming region W75N(B). This region contains two massive protostars, VLA 1 
and VLA 2, separated by <0.7 arc sec [projected separation of <910 astronomical units 
(AU) at the source distance of 1.3 kpc; (P)], both associated with strong water-vapor maser 
emission at 22 GHz ( 10 , 11 ), and with a markedly different outflow geometry. At epoch 
1996, VLA 1 shows an elongated radio continuum emission consistent with a thermal radio 
jet, as well as water maser emission tracing a collimated outflow of ~1300 AU along its 
major axis. In contrast, in VLA 2 the water masers traced a shock-excited shell of ~185 AU 
diameter radially expanding with respect to a central, compact radio continuum source 
(<0.12 arc sec, <160 AU) of unknown nature ( 12 ). A monitoring of the water masers 
toward these two objects from 1996 to 2012 shows that the masers in VLA 1 display a 
persistent linear distribution along the major axis of the radio jet. In the case of VLA 2 
we observe that the water maser shell continues its expansion at ~30 km s -1 16 years 
after its first detection. More important, the shell has evolved from an almost circular 
structure (~185 AU to an elliptical one (~354xl90 AU) oriented northeast-southwest, along 
a direction similar to that of the nearby VLA 1 radio jet ( 13 , 14 ) (fig- SI) and of the ordered 
large-scale (2000 AU) magnetic field observed in the region ( 15 ). The estimated kinematic 
age for the expanding shell is ~25 years, indicating that it is driven by a short-lived, episodic 
outflow event. Moreover, our polarization measurements of the water maser emission show 
that whereas the magnetic field around VLA 1 has not changed over time, the magnetic field 
around VLA 2 changed its orientation following the direction of the major axis of the water 
maser elliptical structure. That is, it now shares a similar northeast-southwest orientation 
with those of both the magnetic field around VLA 1 and the large-scale magnetic field in 
the region ( 14 , 15 )( fig. SI). 

All these observations suggest that we are observing in “real time” the transition from 
an uncollimated outflow to a collimated outflow during the early life of a massive star. This 
scenario predicts that, within the same time span of the evolution of the shell, VLA 2 must 
have also evolved from a compact radio continuum source to an extended elongated source 
along the major axis of the water maser shell. Furthermore, the radio continuum emission of 
VLA 2 should have physical properties (e.g., spectral energy distribution, size of the source 
as a function of frequency) characteristic of free-free emission from a thermal, collimated 
ionized wind ( 16 , 17 ). This can be tested through continuum observations at centimeter 
wavelengths that usually traces the emission from collimated, ionized winds ( 17,18 ). 

Taking advantage of the high sensitivity and high angular resolution of the Jansky Very 
Large Array (VLA) at centimeter wavelengths, we obtained new observations in 2014 at sev¬ 
eral bands in the frequency range from 4 to 48 GHz ( 19 ). These highly sensitive observations 
confirmed the expected scenario proposed above. The source VLA 2 is detected at all bands. 
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In the images of higher frequency bands [U (~15 GHz), K (~23 GHz), and Q (~44 GHz)], 
which have higher angular resolutions (~0.1 to 0.2 arc sec), the source VLA 2 appears clearly 
elongated in the northeast-southwest direction (Fig. 1). Water maser emission was also ob¬ 
served simultaneously with the K band continuum emission (19), allowing a very accurate 
alignment (better than ~1 milli-arc sec) between the masers and the continuum emission. 
We find that the elongation of the continuum emission is in good agreement with that of the 
water maser distribution (Fig. 1). 

Comparison of the radio continuum emission of VLA 2 at K band between epochs 
1996 (10) and 2014 is shown in Fig. 2. Whereas in 1996 the emission was compact, in 
2014 we observed extended emission in the northeast-southwest direction. In particular, 
the core of the radio continuum emission of VLA 2 has evolved from a compact source 
in 1996 (<160 ALT), to an elongated core with a full width at half maximum (FWHM) 
of ~220 AUx<160 AU (position angle ~65°) in 2014. This evolution is consistent with 
those of the expanding water maser shell and of the magnetic field, over the same time 
span (fig. SI) (14)- This core elongation is clearly seen in all the images that we have 
made by weighting the interferometer uv data in different ways, from natural to uniform 
(19) (fig. S2). In addition, we also found that the VLA 2 source exhibits weak extended 
emission at distances of ~390 AU southwest from the peak position, which was not detected 
in 1996. Because the sensitivity attained in our 2014 observations (~10 /iJy beam -1 ) is a 
factor of ~15 better than in 1996 (~150 yuJy beam -1 ), it is plausible that the weak extended 
emission was already present in 1996, but not detected because of the lower sensitivity in 
those observations. That the peak intensity of VLA 2 has significantly decreased from 1996 
to 2014, while total flux densities in both epochs have remained similar (table S3), rather 
suggests that the radio continuum emission from VLA 2 arises from a larger, more elongated 
structure in the 2014 epoch (see model below). In addition, the distribution of the water 
masers follows the observed morphology of the continuum emission at both epochs (Fig. 2), 
indicating that they trace a molecular shell surrounding the radio continuum source. We also 
considered the possibility that the extended radio continuum emission in 2014 is actually 
the signature of a new variable protostar. Nevertheless, this hypothesis is rendered unlikely 
by the spatial coincidence of this new radio continuum emission with VLA 2. 

As a result of the wide bandwidths provided by the VLA, we performed an analysis of 
the behavior of the flux density and of the size of the major axis of VLA 2 with frequency 
(19) (table S3 and fig. S4). This analysis was carried out only at U, K, and Q bands, whose 
high angular resolution allows us to completely isolate the emission of VLA 2 from the other 
sources in the field. We find that the flux density (S v ) at the lowest observed frequencies (LI 
and K bands) is dominated by free-free emission with a spectral index aff ~0.61 (S u oc u a ff). 
Moreover, we also find that the size of the major axis of VLA 2 (0) decreases with frequency 
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as a power law (9 oc A 3 ), with an exponent /3 ~ —0.7. The behavior of both the flux density 
and the size of the major axis of VLA 2 is in good agreement with what is expected from a 
thermal, collimated and ionized wind (16,17). The presence of water maser emission and the 
detection of shock-excited SiO emission highly concentrated towards this source (20) further 
supports a very fast wind interaction region rather than an Hll region. The VLA 2 wind 
could contribute, together with the winds from the other young stellar objects (YSOs) in the 
region, to the large-scale (~2 pc) molecular outflow observed in the W75N(B) star-forming 
region (21). 

At the highest observed frequencies (Q band), the emission deviates by ~30% from the 
expected flux densities of the ionized wind at these frequencies, showing in addition a steeper 
spectral index that can be easily explained by the presence of thermal dust emission. The 
presence of thermal dust emission at frequencies higher than 40 GHz is commonly found in 
protostars and is related to the presence of a dusty circumstellar disk around the protostar 
and/or to a dusty envelope. Our data suggest that the spectral index of the dust component 
in the observed wavelength range is a.d U st — 3 to 4 (fig. S4), a value similar to what is 
commonly found in massive protostars (22,23). 

In summary, our observations reveal that the radio continuum source VLA 2 has evolved 
in only 18 years from being a compact source, to become an elongated source in the same 
direction as that of the water maser shell and the magnetic field. Moreover, the character¬ 
istics of the continuum emission and the distribution of the water masers in VLA 2 are now 
consistent with a thermal, collimated ionized wind, surrounded by a dusty disk or envelope. 
We interpret these results in terms of the evolution of an episodic, short lived, originally 
isotropic ionized wind interacting with a toroidal environmental density stratification. 

The interaction of a wind with the surrounding environment produces a two-shock 
structure that travels away from the star. The outer shock accelerates the environment, 
whereas the inner shock decelerates the wind. For young stars with moderate wind velocities 
(few hundred km s” 1 ), both the inner and outer shocks are, in general, radiative and result 
in the formation of a thin expanding shell bounded by the two shocks. In ( 24 ) we model 
this scenario for the case of a steady isotropic wind with mass loss rate M w and terminal 
velocity V w , surrounded by a torus of dense material with a density distribution of the form 
p(R,6) = p 0 /(l + [ R/R c (9 )] 2 ). R is the distance from the star; R c (9) = i? 0 (l + psin 2 0) 
is the latitude-dependent radial size of the core, with 6 the polar angle measured from 
the symmetry axis of the torus; p 0 is the density of the molecular core; and p is related 
to the density contrast between the equator (9 = 7t/2) and the pole (9 = 0) of the toroidal 
environment. This density distribution implies that for R <C R c the density is nearly uniform, 
that is, independent of direction and radial size, whereas for R ^ R c the density decreases 
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as 1/R 2 , which is the behavior of the density of self-gravitating isothermal spheres at large 
distances. 

Our model (24) shows that at small radii (R <C R c ) the shell grows isotropically 
(with no angular dependence), the initial velocity of the expanding shell is the velocity of 
the wind, and the shell decelerates owing to the incorporation of environmental material 
with zero velocity. The expansion at large radii (R 3> R c ) has two important physical 
characteristics. First, there is a strong angular dependence of the expansion velocity of 
the shell, being much larger along the symmetry axis of the torus than along its equator; 
and second, this expansion velocity asymptotically tends to an angular-dependent constant 
value. Thus, this isotropic wind + torus model can explain the distinctive characteristics of 
our radio continuum and water maser observations (the transition from a nearly spherical 
outflow to a collimated, elongated one over a short period of time, and the expansion of the 
masers at constant velocity), because the shell tends to expand faster along the symmetry 
axis of the surrounding torus (Fig. 3). 

In our model (24 j Fig. 3), we identify the semi-major axis of the water maser shell 
reported at different epochs by ( 14 ) with the evolution of the shell radius along the symmetry 
axis of the surrounding torus (0 = 0), and the semi-minor axis of the water maser shell with 
the shell size along the equator (6 = vr/2). Taking this into account, and assuming that 
the radio continuum source VLA 2 traces the emission of the ionized wind that drives the 
shell, we can find a range of possible physical values to fit the observations. In particular, 
we obtain values for the mass loss rate and terminal wind velocity of M w ~ (4-11) xlO -7 M 0 
yr _1 and V w — 110 to 350 km s _1 , respectively. In addition, for the initial epoch of ejection 
of the wind we found t l ~ 1984 to 1985, with uniform density core radial size R 0 ~ 26 to 
29 AU and particle density n 0 — (4 to 6)xl0' cm -3 . 

All the parameters estimated by our model match those expected for a massive proto¬ 
star embedded in a high-density environment. Furthermore, from the thermal dust emis¬ 
sion at 7 mm (after removal of the free-free contamination) we roughly estimate a total 
(dust+gas) mass of > 0.001 M 0 (19). This mass corresponds to an average particle density 
of >2xl0 7 cm -3 in a region of ~0.2 arc sec (the maximum size of the region where we detect 
7 mm continuum emission; Fig. 1). Within the uncertainties of the estimates (19), we find 
that these values are consistent with the parameters of the toroidal gas structure. 

Our data do not allow us to determine what the morphology of the radio continuum 
emission was like at scales <0.1 arc sec in 1996. In particular, we cannot rule out that the 
source was already slightly elongated with a size smaller than the beam in 1996. However, 
1996 corresponds to an epoch ~10 years after the estimated launch time (1983 to 1984). 
According to our model, with that time interval, some deviation from an originally spherical 



wind should already be expected (Fig. 3). The combination of all our radio continuum results, 
together with other important available data, can provide a coherent scenario for this source. 
Our scenario proposes a change in the collimation of the wind, from a poorly collimated 
outflow (maybe initially spherical) to a highly collimated wind, and this is consistent with 
all the available data: a compact continuum source and water masers poorly collimated in 
1996, an elongated radio continuum source and collimated water masers in 2014. We can 
then affirm that a change in the collimation of the wind has taken place in this time span. 
We also find that the ionized wind, the water maser shell, and the magnetic held in VLA 2 
(at scales of few hundred AU) have all evolved spatially in the same sense, sharing now a 
similar northeast-southwest orientation. This is similar to the orientation of the magnetic 
held in the neighboring radio jet VLA 1 and of the “large-scale” (~2000 AU) magnetic held 
in the region (14,15). Therefore, the magnetic held at these “large scales” might control the 
formation of both objects (VLA 1 and VLA 2). In particular, we think that the magnetic 
held could have favored the formation of the toroidal gas structure around VLA 2, with its 
symmetry axis also oriented northeast-southwest, as in our model. 

The proposed scenario can thus explain the singular evolution of an episodic, short-lived 
(~30 years), originally non-collimated outflow, traced by the water masers, into a collimated 
outflow by its interaction with the ambient medium. Our observations and modeling reveal 
that the collimation of these short-lived outflow events from massive protostars occurs at 
relatively large distances from the central star, at radii between ~30 AU (the radial size 
obtained from our model for the initial expansion of the shell) and ~100 AU (the radius where 
the expanding shell turned out from a round to an elliptical shape; fig. SI) [see also (25,26)}. 
These scales are on the order of those predicted for poorly collimated outflows in magnetized, 
massive collapsing cores during the very early stages (~10 3 to 10 4 years), according to the 
magnetohydrodynamic simulations developed by (27). However, these simulations can only 
produce relatively low-velocity (~10 km s' 1 ) poorly collimated outflows, which seem still 
insufficient to explain the larger wind velocities (>100 km s' 1 ) needed to drive the expanding 
motions of the shock-excited water maser shell in VLA 2. We also considered the scenario 
in which a very young spherical compact HII region embedded in an accretion disk begins 
to expand, following the model in (28). In this model, the HII region evolves to a bipolar 
morphology along the symmetry axis of the accretion disk, producing a wind driven by 
the thermal pressure of the ionized gas. However, the relatively low velocity of that wind 
(<30 km s' 1 ) (26) seems also insufficient to shock excite and drive the expanding water 
maser shell. 

In summary, although episodic, short-lived outflows in massive protostars are probably 
related to episodic increases in the accretion rates, as observed in low-mass star formation 
(29), the origin of poorly collimated outflows with relatively high velocities, as observed 



in VLA 2, deserves further theoretical research. Our observations show the “real time” 
evolution of such an originally poorly collimated outflow into a collimated one. This opens a 
new exciting window of opportunity to study how the basic ingredients of star formation (e.g., 
molecular outflow, ionized wind, magnetic field) evolve over the next few years, providing 
insights that may have important implications for our knowledge of the early stages of high- 
mass star formation. We may be on the brink of describing and modeling in “real time” all 
of these rapid changes. 




Fig. 1.— Radio continuum maps of the core of W75N(B) star-forming region obtained with 
the VLA in 2014 at several bands (C, U, K and Q). Left panel shows the C band (5 GHz) image of 
the region enclosing the sources VLA 1, 2, and 3. Right panels show close-ups to the VLA 2 source at U 
(15 GHz), K (23 GHz) and Q (43 GHz) bands. In the K band panel, we also show the positions of the water 
masers (plus symbols) observed simultaneously to the K band continuum, as well as the methanol masers 
(black dots). As seen in the right panels, VLA 2 appears elongated at all frequencies. This elongation is 
similar to that of the maser distribution. Contours are 4, 8, 16, 32, 64, and 128 times the root mean square 
(RMS) of each map: 30 ^Jy beam -1 (C band), 11 /iJy beam -1 (U band), 10 /rJy beam -1 (K band), and 
12 /iJy beam -1 (Q band). The restoring circular beam of each map (shown in the bottom left corner) is 
0.31", 0.15", 0.12", and 0.07" for C, U, K, and Q band, respectively (0.1" ~ 130 AU at the source distance 
of 1.3 kpc). In all panels, the absolute position of the (0,0) is at right ascension RA(J2000) = 20 ,t 38 m 36.486 s 
and declination DEC(J2000) = +42°37'34.09" (± 0.03"), the peak position of VLA 2 at Q band, where the 
massive protostar is expected to be located. The accuracy in the relative positions of the water and methanol 
masers with respect to the K band continuum is better than ~1 and ~30 milli-arc sec, respectively ( 19 ). 
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Fig. 2.— Comparison of the K band continuum emission of VLA 2 in epochs 1996 (top) and 
2014 (bottom). The (0,0) position is the same as in Fig. 1. Contours are 30, 40, 50, 60, 70, 80, and 90% of 
the peak intensity in 1996 (1.42 mJy beam -1 ) (RMS = 150 /rJy beam -1 ), and 5, 10, 20, 30, 40, 50, 60, 70, 
80, and 90% of the peak intensity in 2014 (0.82 mJy beam -1 ) (RMS = 10 /iJy beam -1 ). Both maps were 
obtained with the same restoring circular beam of 0.12" (shown in the bottom left corner of each panel). 
In both panels, the half-power level is shown as a dashed line. We also show the water maser positions 
(plus symbols) for epochs 1996 and 2014 as observed with the VLA by ( 10 ) and this work, respectively. 
The position of the methanol masers for epoch 2014 (black dots; this work) are also indicated. The FWHM 
size of the radio continuum emission has evolved from a compact source into an elongated source in the 
northeast-southwest direction, in a direction similar to that of the observed evolution of the water maser 
shell and magnetic field in VLA 2 ( 14 ) (fig- SI). 
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Fig. 3.— Evolution of an isotropic ionized wind driving the shock-excited water maser shell 
and that interacts with a toroidal environmental density stratification ( 24 )• The model predicts 
an initial round shell, evolving into an elongated shell that expands faster along the symmetry axis of the 
torus. Here we represent density isocontours of the toroidal environment (blue contours) for values p/po = 
0.1, 0.3, 0.5, 0.7, and 0.9 [eqs. El and E2 with p = 5; (24)]- The spatial evolution of the shell as a function 
of time t after the initial ejection of the wind (red contours) is shown for values t/to = 1, 5, 10, 20, 30, and 
40 with to = 0.5 years [adopting V w = 250 km s -1 and Rq = 28 AU from (24 )]• 
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